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Abstract 

The 5D Kepler system possesses many interesting properties. This system is superintegrable and 
also with a su{2) nonAbelian monopole interaction (Yang-Coulomb monopole). This system is also 
related to a 8D isotropic harmonic oscillator by a Hurwitz transformation. We introduce a new 
superintegrable Hamiltonian that consists in a 5D Kepler system with new terms of Smorodinsky- 
Winternitz type. We obtain the integrals of motion of this systems. They generate a quadratic 
algebra with structure constants involving the Casimir operator of a so(4) Lie algebra. We also 
show that this system remains superintegrable with a su{2) nonAbelian monopole (generalized 
Yang-Coulomb monopole). We study this system using parabolic coordinates and obtain from 
Hurwitz transformation its dual that is a 8D singular oscillator. This 8D singular oscillator is also 
a new superintegrable system and multiseparable. We obtained its quadratic algebra that involves 
two Casimir operators of so(4) Lie algebras. This correspondence is used to obtain algebraically 
the energy spectrum of the generalized Yang Coulomb monopole. 
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I. INTRODUCTION 



In classical mechanics a Hamiltonian system with Hamiltonian H and integrals of motion 



is called completely integrable (or Liouville integrable) if it allows n integrals of motion 
(including the Hamiltonian) that are well defined functions on phase space, are in involution 
{H,Xa}p — 0, {Xa,Xh}p — 0, a,b — l,...,n — 1 and are functionally independent {{,}p is a 
Poisson bracket). A system is superintegrable if it is integrable and allows further integrals 
of motion Yi){x,p), {H,Yi)}p = 0, b=n,n+l,...,n+k that are also well defined functions on 
phase space and the integrals{//, Xi, Yn, Yn+k} are functionally independent. A 
system is maximally superintegrable if the set contains 2n — 1 functions and minimally 
superintegrable if it contains n + 1 such integrals. The integrals Yj, are not required to be 
in involution with Xi,...Xn-i, nor with each other. The same definitions (i.e. integrability 
and superintegrability) apply in the context of quantum mechanics but {H, Xa, Yi,} are well 
defined quantum mechanical operators, assumed to form an algebraically independent set. 
The most well known superintegrable system is the Kepler-Coulomb potential [1-3]. Its 
integrals of motion consist of the angular momentum and the Laplace-Runge-Lenz vector 
whose origin can be traced back 80 years earher [4-6] than Treatise of celestial mechanics of 
Laplace [7] with the works of BernouUi [6] and Hermann [4,5]. 

A systematic classification of classical and quantum superintegrable systems was initiated 
40 years ago by Smorodinsky, Winternitz and co-workers for Hamiltonians with quadratic 
(in momenta) integrals of motion in two-dimensional Euclidean space [8]. Four types of 
maximally superintegrable systems were obtained with multiseparability [8]. This search 
was extended in three-dimensional Euclidean space, where five potentials are maximally 
superintegrable with five integrals of motion. Nine potentials belong to the class of minimally 
superintegrable systems and are characterized by the existence of four second order integrals 
of motion [9-11]. One of these nine systems include the generahzed Kepler system [12-23] 



X, 



a 




a = 1, n — 1 





2 r r(r -|- xq) r{r — xq) ' 



where r = •\/xf~+xf+~x^ and ci and C2 are non-negative constants. When c\ — ci ^ 0, 
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this system reduces to the Hartmann potential introduced in 1972 as a model for ring- 
shaped molecules [24] . As the superintegrabihty property is preserved for the Kepler system 
with a monopole interaction (MICZ-Kepler) [25-31], the generalized Kepler system with 
monopole (generalized MICZ-Kcpler system) is also superintcgrablc [32-36]. This system 
is also related to another superintegrable systems, a 4D singular osicllator [33,36-38], by 
a Kustaanheimo-Stiefel transformation [39]. Let us mention that the 2D Kepler system is 
related to the 2D isotropic harmonic oscillator by a Levi-Civita transformation [40]. This 
relation that map these two superintegrable systems in classical and quantum mechanics 
was introduced at the origin in context of the regularization of the classical Kepler problem. 
A such relation also exists for the 5D Kepler system that can be related to a 8D isotropic 
harmonic oscillator by Hurwitz transformations [41]. This system and the Yang-Coulomb 
monopole ( i.e. 5D Kepler system with a su{2) non-Abelian monopole) have attracted 
a lot of attention [42-51]. Let us mention that the existence of the Levi-Civita and its 
generalizations is one important property of these systems. The Hurwitz relations were 
exploited to present a derivation of the energy spectrum of the Yang-Coulomb monopole 
from the ladder operators of the 8D harmonic oscillator [50]. There is an analog of these 
relations between ID harmonic oscillator and Kepler system. These four cases are related 
to the Euler identity (xq + ... + x'^_-^ = {u^ -|- ... -|- ujy_i)^) which has a bilinear solution in Ui 
only for the following pairs (D,d): (1,1), (2,2), (4,3) and (8,5) [47]. 

Quadratic algebras and more generally polynomial algebras have been applied to many 
quantum systems [22,23,36,52-78]. Such algebraic structure generated by the integrals of 
motion do not only provide an elegant method to obtain the energy spectrum but also 
explain the degeneracies of the energy spectrum. Recently [36,70,72], it was pointed out 
how results obtained in context of 2D superintegrable systems and their quadratic algebras 
can be applied to three and four-dimensional superintegrable systems. 

The purpose of this paper is to study a generalized 5D Kepler system and moreover to 
obtain the superintegrable analog with a nonAbelian monopole. The study of systems with 
nonAbelian monopole with Smorodinsky-Winternitz terms and the related duality trans- 
formation are unexplored subjects. The paper is also devoted to the study these systems 
using representation theory of quadratic algebra developped in context of 2D superintegrable 
systems. 

Let us present the organization of the paper. In Section 2, we recall results obtained in 
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context of 2D systems concerning the realizations of quadratic algebras as deformed oscillator 
algebras and the representation theory. In Section 3, we recall results concerning the five- 
dimensional Kepler system and introduce a generalized 5D system. We obtain its integrals 
of motion and we show that they generate a quadratic algebra with structure constants 
involving an integral that commute with the generators of this quadratic algebra. This 
integral is also the Casimir of a so(4) Lie algebra. We present an algebraic derivation of 
the energy spectrum of this system. In Section 4, we introduce a generalized Yang-Coulomb 
monopole, present its integrals of motion, the energy spectrum and the wavefunctions using 
parabohc coordinates. We will show also that this system is related to a 8D singular oscillator 
by a Hurwitz transformation. We will study the quadratic algebra of this system and present 
the finite dimensional unitary representations. This quadratic algebra involve two integrals 
that commute with the generators that are the Casimir of so(4) Lie algebras. We also use 
these results to reobtain the energy spectrum of the generalized Yang-Coulomb monopole 
algebraically. 

II. QUADRATIC ALGEBRAS 

Quadratic and more generally polynomial algebras were introduced in the context of non- 
relativistic quantum mechanics [22,23,55-59] 20 years ago. Their finite-dimensional unitary 
representations can be used to calculate algebraically the energy spectrum of superintegrable 
systems. Let us recall results obtained in context of two-dimensional quantum superinte- 
grable systems. We refer the reader to the Ref.60 for more detailed discussions. 

In the case of a two-dimensional Hamiltonian H allowing two second order integrals A 
and B (i.e. [i/, A\ = [H, B] = 0), we have the following graph ( as in Ref. 72 ) 



The most general quadratic algebra generated by these integrals is given by the following 
commutation relations [60] : 



A 



H 



B . 



(3) 



[A, B] = C, [A, C] = aA^ 



+ 7{A, B} + 5A + eB + (, 



(4) 



[B, C] = a^' 



- a{A, B} + dA-SB + z. 
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The structure constants 7, S, a, e, ^ and z of the quadratic algebra given by Eq.(4) are 
polynomials of the Hamiltonian H. The degree of these polynomials can be deduced from 
the order of these commutators. The structure constants a, 7 and a are constant and 



S = Sq + SiH, e = eo + eiH, d = do + diH, (5) 

C = Co + Cl^ + C2^^', Z^Zo + ZiH + Z2H^ 

The Casimir operator (i.e. [K, A] = [K, B] = [K, C] = 0) of this quadratic algebra is 
thus given in terms of the generators {A,B and C) by 

K = C''- a{A\ B} - ^{A, B^} + {a^ - 6){A, B} + (7^ - e)B^ (6) 

+ (75 - 2C)B + ^A' + {d+^ + a')A^ + (^ + aS + 2z)A. 

In order to obtain the finite-dimensional unitary representations, we consider realizations 
of the quadratic algebras in terms of deformed oscillators algebras [56] {N,b'',b} satisfying 
the following equations : 

[TV, 6t] ^ {,t^ [jv, b] = -b, bb^ = + 1), b^b = ^N), (7) 

where ^{x) is a real function called the structure function satisfying $(0) = and > 
for a; > 0. We have the existence of finite dimensional unitary representations when we 
impose the existence of a positive integer p such that + 1) = 0. The reahzations of the 
quadratic algebra given by Eq.(4) of the form {A = A{N), B = b{N) + b'^p{N) + p{N)b) were 
studied [60]. Two cases were obtained (i.e. 7 7^ and 7 = 0). Let us only present the results 
in the case 7 7^ with the following choice of structure constants {a — a — 5 — 0) similarly 
to the generahzed MICZ-Kepler systems ( it is a constant determined from the constraints 
on the structure function ): 

^^^^ " 2i23^8(Ar + u){l + N + u){l + 2{N + u))'' 

$(Ar) = -3072-f^K(-l + 2(N + u))^ + A8d'y^(-3 + 2(N + u))(-l + 2{N + u))^(l + 2{N + u)) 
+I22887V + 327^(-l + 2{N + u)y{-l - 12{N + u) + 12{N + uy){-4:d€-f^ + 87^^) 
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-2567^(-l + 2(N + M))^(-3de^7^ + 2rfe7^ + Ue^^z - A-f^z). 

The Casimir operator K can be written in terms of the Hamiltonian only. We have a 
energy dependent Fock space of dimension p+1 if 

^p+l,u,E)^0, $(0,M,E) = 0, $(x)>0, Vx>0 . (9) 

The Fock space is defined by 

H\E,n>= E\E,n>, N\E,n >= n\E,n > b\E,0>^0. 

b^\n >= ^/^n + l,E)\E,n + l >, b\n >= y/<5{n, E)\E,n - 1 > . 

The energy E and the constant u are solutions of the equations obtained by the system given 
by Eq.(9). They represent the finite-dimensional unitary representations with dimension 
p+1. 

A. Higher dimensional systems and quadratic algebras 

These results obtained in context of two-dimensional systems can be applied to various 
systems in three and four dimensions. The Kepler system [55], the harmonic oscillator 
[56] and the MICZ-Kcplcr systems on three sphere [73,31], the Hartman potential [23], the 
generalized Kepler [22], generalized MICZ-Kcplcr systems [36], the 4D singular oscillator 
[36], the generalized Kaluza-Klein monopole [70] and the non-degenerate extended Kepler- 
Coulomb system [72] possess a quadratic algebra with three generators. 

We allow in context of higher dimensional superintegrable system the structure constants 
a, S, 7, e, (, a, d and z to be polynomial functions not only of the Hamiltonian but also of 
any other integrals of motion (Fj) that commute with the Hamiltonian, with each other and 
also with the generators of the quadratic algebra A, B and C (that is, [Fi,H] = [Fi,Fj] = 
[F,,A]^[F,,B]^[F,,C]^Q). 
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A B Fi Fj (10) 

When we study this algebra's realization in terms of deformed oscillator algebras and 
its representations, we will fixed the energy (Hip — Eip) but also these other integrals 



{Fi^l; — fill)) of motion that form, with the Hamiltonian, an Abehan subalgebra. The 
Casimir operator K of this quadratic algebra is thus given in terms of the generators and 
will be also rewritten as a polynomial of H and F^. 



III. GENERALIZED FIVE-DIMENSIONAL KEPLER 

A. Five-dimensional Kepler system 

Let us first recall some known results concerning the 5D-Kepler system given by : 

2 r' ^ ' 

with r = ^Xq + x\ + x\ + x\ + xl and 

P^ = -ihdi, [xi, Pj] = iMij. (12) 
This system has the following integrals of motion ([if, L^j] — [H, M^] — 0) 

Lij ^ XiPj - XjPi, Mk^^{PiLik + LikPi) + ^^. (13) 
These integrals generate the following dynamical symmetry algebra : 

[Z/jj, Z/jyjji] ihSiraLjfi ifi6jimLiii ihdinLjul -\- ihSjnLiYm (-^'^) 

[Lij, Mk] = ihSikMj - ihSjkMi, [Mi, Mk] = -2ihHLik. (15) 

The relations given by Eq.(14) generate a so{5) Lie algebra. The algebra generated when 
we consider also the Eq.(15) and we fixed the energy is isomorphic to the so(6) Lie algebra 
for the bound states. The symmetry algebra can also be seen as a Kac-Moody algebra [79] . 
This dynamical symmetry algebra can be used to obtain algebraically the energy spectrum. 

B. Generalized 5D Kepler system 

Let us now introduce the following Hamiltonian : 

if = + + (16) 

2 r r(r + xq) r{r — xq) 
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The two last terms break the so(5) rotational invariance. This systems is superintegrable 
and has the following integrals of motion : 

^ = iV^ + ^, s = M„ + £f^-£f±^, (17) 

(r + Xo) [r — xq) r[r + Xq) r{r — Xq) 

— -^ij) ^ij — ^iPj ~ ^jPh (18) 

l<i<j<4 

where i,j=l,2,3,4. 

We have the commutation relations : 

[H, A] = [H, B] = [H, L'] = [A, L'] = [B, L'] = [H, L,,] = 0, (19) 

and in terms of graph : 



/ 

/ 

/ 

/ 

/ 

/ 

/ 



H H 

\ / \ 

\ / \ 

\ / ^ 



A B (20) 

These integrals generate the following quadratic algebra 

[A,B] = C, [A,C] = 2h^{A,B} + Sti^B-A{ci-C2)h^co, (21) 

[S, C] = -2h^B'^ + ^h^HA - Ah^L^H + IQh^H - 8h^{ci + C2)H + 2h^cl, (22) 

K = im^HL"^ - 8h^(ci - C2fH + 32(ci + C2)h'^H - 32ffH (23) 
+4h^clL'^ + 8h^{ci + C2)cl - Ah^cl. 

This quadratic algebra has a very interesting structure. The structure constants involve 

the operator that is the Casimir operator of the so(4) Lie algebra generated by Ljj for 
i,j = 1,2,3,4. The / ( L'^ip = lip) can be calculated from the representation theory and 
Casimir operators of so(4) Lie algebras. 

C. Finite dimensional unitary representations 

We obtain for the generalized 5D Kepler system the following structure function : 
$(x) = 619U56Eh^^{x+u-{^{l-mi-m2))){x+u-{^{l-mi+m2))){x+u-{^{l+mi-m2))) 



(x + u- (^(1 + mi + m2))){x + u-{^ - ^-^^ 



Co xw , A 

))[x + u 



Co 



:)) (24) 



2 ^/^2Eh' 

with ti?m\ 2 = 16ci,2 + 4/ + 4fi^. From the constraints on the structure function given by the 
equation (9) we obtain u—\-\- ^^^^^ and 



;i2(p + i + m±m)2^ 



p = 0, 1,2,3,... 



(25) 



, 6291456n^Sx(p+ 1 - x)c^, , . ^^ ^ /^.x 

$ X = , ^ ■ f^(p + l + mi-x)(p+l + m2-a;)(p + l + mi + m2). (26) 

+ 1 + mi + 7712) 

IV. GENERALIZED YANG-COULOMB MONOPOLE SYSTEM AND DUALITY 
A. Yang-Coulomb monopole 



The Yang-Coulomb monopole is given by the following Hamiltonian that consists in a 
5D-Kepler systems with a su{2) monopole interaction : 

with TTj = —ihdxj — hA'jTa and 



r{r + xO) ^ 



(27) 



rr, / / \ f n \ ^ / \ ^ COsiaT) d 

actT c/pT sin{pT) ajT 
±2 = t{stn[aT)cot[pT) ^ cos(aT)^-5 . . ^ \ "a — /' 



sin{/3T) d^x' 



d 



daT ' 



(28) 

(29) 
(30) 



with (i,j=0, 1,2,3,4 ; //,z^=l,2,3,4).These operators given by Eq.(27)-(30) satisfy the following 
relations : 



(31) 



We have : 


















—ia 





ia^ 






'^O ^ 



r — - 
2 



/ 



ia^ 
-ia^ 



^0 0^ 



a2 

(T^ 



(32) 



where cr* for i = 1, 2, 3 are the Pauli matrices : 



a 



1 

1 



a 



-i 

1 



a 



1 
-1 



(33) 
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We can construct the following field tensor 



Ft,^d,At-dkAt + eai,cA'lAl (34) 
and using Eq(28)-(30) and (34) generate the following integrals of motion : 

Lik^iXiTTk-XklTij-r'^hF^lTa, Mk ^ -{TTiLik + LikTTi + 2co—) (35) 

They satisfy Eq.(14) and (15) and thus preserve the so(6) Lie algebra. 

B. Yang-Coulomb monopole with Smorodinsky-Winternitz terms 

We can introduced the following Hamiltonian : 

H = -(TT,)^ + - ^ + + (36) 

2 2r^ r r[r + xq) r[r — xq) 

Similarly, we can show with more involving calculations that this system allows integrals 
of the form given by equation (17) and (18) where the Lij and Mq are given by Eq.(35). 
This system is thus a new superintegrable system and multiseparable ( a consequence of 
this existence of more than one second order integrals of motion). It allows separation of 
variables in hypcrsphcrical, spheroidal and parabolic coordinates systems. Let us present the 
wavefunctions and a derivation of the energy spectrum of this systems using the parabolic 
coordinates. We define the parabolic coordinates by : 

= ^(/w-i^), (37) 

. »(«+7) 

Xi + 1X2 ^ y/Jjil'COS{-)e 2 , 

• f P \ '("~t) 

x^ + 1X4^ = ^ jiusiny — je 2 ^ 

with ly, & [0, 00). The Schroedinger equation given given by the Hamiltonian (36) (with 
Hip — eip) become (with Ji — Li + Tj) 

(4j2-2ci) (4L2-2C2),, 2, 2co , , ^ 
fj,{li + v) -\-v) fJ, + v 

with 
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We use the following form for the wavefunctions (similarly to the Yang-Coulomb monopole 
[44]) 

i; = Ch{fx)Mu)Di^^,,,{a,/3,-f,aTjT,lT), (40) 
with Di^^,^,{a, /3, 7, ar, Pt, It) given by 



Di^m'Ac^, P, 7, «T, /3t, it) = ^ (2-^+lK2r+l) ; ^ C/f^^^,D^^,(a, ^, 7) AT'(«t, ^t, 7t), 

M=m+t 

(41) 

where C'l^.rp^ are the Clebseh-Gordon coefficients, D^^,{a, j3, 7) and Dj^,{aT, Pt, 7t) are 
the Wigner functions, and C is a normalization constant. Using the following parameters 
si — 4(j(j + 1) — 2ci), S2 — 4(/(Z + 1) — 2C2), (5 = |l and a = ^ and the wavefunctions as 
given by Eq.(40), we obtain two seconder order differential equations : 

{d,{f,d,) -^ + ^ + ^f^)fM = IfM, (42) 

{d.iud^) - £ + f + f ^)/2 W = -^/2(^). (43) 

The solution can be written in term of the hypergeometric function (with x — n^—^ 

and y = uy^—f^) 

fi{n) ^ x^e~^F{-ni,si + l,x), fiiiy) ^ y^e'^ F{-n2, S2 + l,y), (44) 
where ni, n2 — 0, 1, 2, ... and 

si + 1 (2w-a) S2 + I {2v + a) , 

ni = h ^ — , ^2 = ^ — , (45 

2 4v^ 2 4v^ ^ ^ 

From Eq.(45) we have the relation 

ni + n2 + (si + S2 + l) = -7-^, (46) 
that allows us to obtain the energy spectrum 



-Co 



2h^{ni + n2 + {si + S2 + I))'' 



(47) 
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C. Hurwitz transformation and 8D singular oscillator 



The Schroedinger equation of the generahzed Yang-Coulomb monopole system 
is related to the following 8D singular oscillator 

by the following Hurwitz transformation (that convert to the direct product ® S^) 



xo = + + M3 - - M5, (50) 

Xi = 2(mo'U4 — — 112^6 1^31*7)) 

X2 = 2(woM5 + 1*1^4 - M2l*7 + 1*31*6), 
X3 = 2(Moti6 + 1*11*7 + 1*2l*4 - W3M5), 
X/^ = 2(^0 1*7 - 1fl1*6 + 1*21*5 + 1*31*4), 

1 («o + mi)(^2 - W3) ^ 

2 (Uo - ZMl)(l*2 + i1*3) 

^2 _|_ ^2 2 

Bt — 1 arctan — ^ \ , e [0, tt] 

1^2 + ^3 

^ (Mo-«Mi)(i*2-i«3) ^ rn . N 

with the following relations between the parameters 

E —up" 

Co = J, e = Cj = 2Ai. (51) 

The 8D singular oscillator is also superintegrable and has the following integrals of motion 



Jii = Uidu^ - Ujdui , i<j, z, j = 0, 1, 2, 3 (52) 

Kij = Uiduj - Ujd^^, i < 3, i, j = 4, 5, 6, 7 (53) 

E 4' E 4' (54) 

l<i<j<4 l<i<j<4 
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9^ 2 , 2 , 2 , 2 2 2 2 2\ , -^1 ^ 



2 ^9^2 ' 2 ^"'^ ' ' ' ' „2 + ^,2 + „2 ^ ^,2 ^2^„2 + „2 + ^2- 

(56) 

We have the following commutation relations : 

[H, A] = [H, B] = [H, K'] = [H, j'] = [H, J,,] = [H, K,,] = (57) 

[A, J'] = [A, K'] = [B, J'] = [B, K'] = 
and in terms of graph : 

H H 

//\\ / \ 

/ \ \ ^ ^ 



/ / ^ \ 

/ , \ i^2 K,, 

A I \ B 

\ "x/ )^ I H 

\ ^ / > / / \ 

\ / X 1 / / N 

\ / /^x * / 

\ / / w 

\ / / \ I / / \ 

j2 ^2 J2 J.. (58) 

The integrals generate the following quadratic algebra 

U, B] = C, [A, C] = 2{A, B} + 8B + J^H - K^H - ^i^i-lM//^ (59) 

nr 



[5, C\ = m^B'^ + 2i7^ - im^uj'^A - 4hWj^ - AHWk'^ + 8(Ai + A2 - 4n^)a;^ (60) 



K = -2 J^i/^ _ 2K^H^ + 4(Ai + A2 ^V ^^ ^ ;^2^2 ^ ^2^2 _ ;^2^2(^2)2 _ 2h^u\PK^ (61) 



-4(A. - A. - AteWJ^ + 4(A, - A. + 4n>^X^ + 4((A. - A.)- - 8(A. + A.)/.- + 16;^V_ 

The structure constants and the Casimir operator allow us to obtain the structure function 
$(x) of the deformed oscillator algebra. The energy and the two Casimir operators are 
fixed, i^V — Eil^: J'^fp — ji^ ^iid -^^V^ — kip. The structure function is then given by (with 
h'^rrii = h^f + 2Ai + /i^ and /i^ms = h'^k'^ + 2A2 + h"^) 
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^ 3 . 2-.= (. + „ _ ( 1 _ ) (. , . - ( 1 + I , (62) 

x + u- Q - ^(-mi + ma)^^ + m - Q - ^(-"^1 - ms) 

We thus found the following energy spectrum with the corresponding finite dimensional 
unitary representation: 

4>{x) = 3 • 2^^uj'^x{p + 1 — x){p + 1 + nil — x){p + 1 + 1712 — x) (64) 

(p + 1 + mi + 7712 - x){2p + 2 + mi + m2 — x). 

Using the relations between the parameter of the generalized Yang- Coulomb monopole 
and the 8D singular oscillator, we obtain the following equation (with Eq(51)) 

4co = E = 2u;hip + 1 + !I!l±I!^) ^ 2^eh{p + 1 + ^^^4^), (65) 

that can be used to reobtain algebraically the energy spectrum of the generalized Yang- 
Coulomb monopole 

This results coincide with the one obtained from a direct approach using the parabolic 
coordinates and given by Eq.(47) ( with m^ = 2sj and p — ni + n2 )■ 



V. CONCLUSION 



One main results of this paper is the study of three new superintegrable systems the 
generalized 5D Kepler system, the generalized Yang-Coulomb monopole and the 8D singular 
oscillator. The study of systems with nonAbelian monopole is important, and this paper 
point out how they could be generalized using Smorodinsky-Winternitz terms. We show 
that for the generalized 5D Kepler systems (as for the well studied 5D Kepler systems) the 
superintegrability property is preserved when one add the appropriate nonAbelian monopole. 
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The Hurwitz transformation maps also this superintegrable system to an other one a 8D 
singular oscillator. 

The algebraic structure presented in this paper are interesting. In the case of the gen- 
eralized Kepler system the quadratic algebra involve the Casimir operator of a so(4) Lie 
algebra. In the case of the dual of the generalized Yang-Coulomb monopole, the 8D singular 
oscillator, the quadratic algebra involve two integrals that are Casimir operator of so(4) Lie 
algebras. This point out also how the quadratic algebra approach can be apply in context 
of higher dimensional superintegrable. 

Let us mention that these results can have a wider applications, it was shown that sys- 
tems with monopole can be extended in N dimensions [80] and a generalized Hurwitz trans- 
formation were obtained for a 9-dimensional hydrogenlike system with so(8) non-Abelian 
monopole [81]. Let us mention also recent works on monopoles [82-86]. 
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